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Abstract. The role of geometry on dispersive forces is investigated by calculating the 
energy between different spheroidal particles and planar surfaces, both with arbitrary 
dielectric properties. The energy is obtained in the non-retarded limit using a spectral 
representation formalism and calculating the interaction between the surface plasmons 
of the two macroscopic bodies. The energy is a power-law function of the separation 
of the bodies, where the exponent value depends on the geometrical parameters of the 
system, like the separation distance between bodies, and the aspect ratio among minor 
and major axes of the spheroid. 



1. Introduction 

The Casimir effect is one of the macroscopic manifestations of the fluctuations of the 
quantum vacuum pQ. Casimir showed that the energy K(z), between two parallel perfect 
conductor plates can be found from the change of the zero-point energy of the classical 
electromagnetic field, as 



where Ui(z) are the proper modes that satisfy the boundary conditions of the 
electromagnetic field at the plates which are separated a distance z. The energy obtained 
by Casimir is a power law function of z, and at large distances U{z) oc z~~ 3 , while at 
short distances U{z) oc z~ 2 . Later, Lifshitz obtained a formula to calculate the force 
between two parallel half-spaces with arbitrary dielectric properties [2]. The Lifshitz 
formula depends only on the reflection amplitude coefficients of the half-spaces and 
the separation between them, finding the same dependence of the energy with z. In 
1968, van Kampen et al. |3] showed that the Lifshitz formula, in the non-retarded limit, 
is obtained from the zero-point energy resulting from the Coulomb interaction of the 
surface plasmons of the plates. After, Gerlach J3j did an extension showing that also in 
the retarded limit, the Lifshitz formula is obtained from the interacting surface plasmons. 
Recently, it has been shown that the Casimir energy is given by the contribution from 
the interacting surface plasmons and propagating modes in the cavity formed by the 
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parallel plates jHj- Furthermore, it has been shown that the contribution from the surface 
plasmons is essential to calculate the Casimir energy [Sj. 

Surface plasmons are evanescent electromagnetic waves that propagate along the 
surface of conductors, and vanish elsewhere. By altering the surface, i.e., modifying 
size, shape, and/or environment of the conductor, the properties of surface plasmons 
can be tailored jOJEl- Therefore, if the shape of at least one of the bodies is modified, we 
would expect to observe changes on the energy due to the interaction between surface 
plasmons. In this paper, we study the influence of the geometry on the zero-point energy 
due to the interaction between macroscopic bodies. Using a method based on a spectral 
representation formalism [E], which determines the proper frequencies of the system, 
we calculate the zero-point energy of all the interacting surface plasmons between a 
spheroidal particle and a flat plate, both with arbitrary dielectric functions. We find 
that the geometry plays an important role in the determination of the energy. 

2. Interacting surface plasmons 

We consider a spheroidal particle located near a flat substrate. The particle is generated 
by the rotation around one of the axes of an ellipse with lengths 2r> and 2r<, with 
r > > r<. The symmetry axis of the particle is perpendicular to the substrate, and its 
center is located at a distance d from the substrate which has a dielectric constant e su b- 
The particle has a dielectric function e part and is embedded in an ambient of dielectric 
constant e am b, which is equal to 1 for vacuum. We consider that the three media: 
particle, substrate and ambient are non-magnetic. We are not considering nonlocal 
effects, therefore, all the variables are function of the frequency only. The explicit 
dependence of the following equations with the frequency is omitted here for simplicity. 
The charge distribution on the particle's surface, in the presence of a substrate, depends 
with the components of the external electromagnetic excitations, throughout the so 
called effective polarizability tensor a e R. When the particle is far from the substrate 
the effective polarizability becomes the polarizability of the isolated particle. But, when 
the particle is close to the substrate the induced multipolar interactions modify the 
electromagnetic response of the system because more surface plasmon resonances of the 
particle interact with those in the substrate [T5j . 

The quantum vacuum fluctuations induce a charge distribution on the particle 
which also induces a charge distribution in the substrate. In the non-retarded limit, the 
induced lm-th multipolar moment on the particle is given by [Hj 



where V^ c (u) is the field associated to the quantum vacuum fluctuations at the 
zero-point energy, V^^u) is the induced field due to the presence of the half-space, 
and a l ^(u) is the lm-th component of the effective polarizability of the particle. It 
is known that the poles of a e g yield the frequencies of the proper modes of the 
system [H1[TUJ[TT]. In the non-retarded limit, the spheroidal radius of the major axis, 
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and the minimun separation between the particle and the substrate are much smaller 
than the characteristic length of the system given, in this case, by I = c/uj p , with c the 
speed of light, and u> p the plasma frequency of the metallic particle. 

The analysis of a e fj for the system described above was done as follows. First, 
the electric potential induced in the system at any point in space was calculated 
to all multipolar orders. To find the solution for the induced potential a spectral 
representation of the Fuchs-Bergman-Milton type [HJ EH E] was developed [12]. By 
identifying the multipolar moments Qimiu) induced in the particle, the components of 
a e s were obtained. Then, the frequencies of the proper modes for different shapes and 
locations of the particles is calculated by choosing a model for the dielectric function of 
the particle. For a detailed description of the method see Ref . [T3] . 

Within the spectral representation formalism, we can write the component of a e g 
in the following form: 



a 



Im i 
eff ' 



s,q 



sq 



where v is the volume of the particle, z is the minimum separation distance from the 
particle to the substrate; such that, z = d — r> for prolate spheroids while z = d — r< 
for oblate spheroids. Here, 

U(u) = [1 - epart/eamb]" 1 , (4) 

is the so-called spectral variable and the strengths C l ™ = {U 1 ™) 2 are the so-called 
spectral functions where U ™ is the unitary matrix that satisfies the relation 

(U^Y l Htt(z)U%=4™%{z). (5) 

The matrix Hf,^,^) depends only on the geometrical properties of the model and on 
the dielectric properties of substrate and ambient, through the contrast parameter 

fc (eamb ^sub)/(^amb ^sub )• (6) 

Furthermore, H*,^,(z) is a real and symmetric matrix which is given by [H]: 

HZ(z) = <{z - oo)^ m + f c D l ™(z), (7) 

here n l ™(z — > oo) are the depolarization factors of an isolated spheroid, and D l ™(z) is 
a matrix given by the multipolar coupling due to the presence of the substrate; this 
later vanishes when z — > oo. Note that H^^z) contains all the information of the 
geometry of the system and the dielectric constant of the substrate, and is independent 
of the dielectric properties of the particle. In Eq. ©, a e s(uj,z) is given as the sum of 
terms which show resonances at frequencies uj given by the poles of the equation, i.e., 
when u{u) = n l ™(z). As a consequence, an explicit procedure to calculate the strength 
and position of the resonances was obtained. In the next section we present results for 
the energy between a spheroid and a plate due to the interaction of all the induced 
multipolar surface plasmons. 
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3. Zero-point energy of a spheroid near a substrate 

As follows, we present a systematic study of the energy for oblate and prolate particles 
with different asymmetries and considering different kinds of substrates. The energy is 
calculated by substituting the frequencies of the proper modes, obtained as described 
above, in Eq. (JTJ). Using the spectral representation, a systematic study of the energy 
in term of the geometrical parameters: r>, r<, and z, and the dielectric properties of 
the substrate, can be done. 

In the presence of a dielectric half-space, the proper modes, given by a e s(uj, z), are 
red-shifted always as the particle approaches the substrate, and this shift depends on 
the separation distance z. In general, the interaction energy is negative for any z, and is 
proportional to (1 + z)' 13 , with f3 = 2L + 1. Here L is a positive integer which labels the 
highest order of the multipolar interaction, and depends on the geometrical parameters 
r>, r<, and z fH]. For example, when r > = r< = a, and z > 5a, then the relevant 
multipolar excitations are given by L = 1, i.e., only surface plasmons with dipolar 
distributions are important. On the other hand, if 5a > z > 2a the interactions among 
surface plasmons with dipolar and quadrupolar charge distributions become relevant 
and L = 2, such that, the energy, when z ~ 2a, is proportional to (1 + z)~ 4 . As z — > 0, 
more and more multipolar charge distributions must be taken into account, and when 
the spheroid is touching the substrate, one would expect that L — > oo, so that also the 
interaction energy also does. 

As study, we employ the plasma model for the dielectric function of the 

particle, e part (cG>) = 1 — where u p is the plasma frequency, which is different for 
different metals. Therefore, the frequencies of the proper modes, obtained from the 
relation u{uf) = n l ™(z), are 



u%(z) =cj, jn%(z), (8) 
and according to Eq. (JTJ, the zero-point energy is given by 

U{z) = [yfW - ^<(* - oo)j , (9) 

sq,lm 

where n l ™(z — > oo) denotes the proper modes of the isolated spheroid. From eq. El we 
can define a dimensionless energy H = U(z)/fkJ p , and study in detail the behavior of 
the system independently of the plasma frequency of the metallic spheroid. 

Let us first analyze the influence of the dielectric properties of the substrate on 
the energy as a function of the dimensionless distance z/r < . In Fig. Q we show 
— log(— H) for an oblate spheroid with an aspect ratio among the semi-axes of the 
spheroid, r > /r < = 1.4, and different substrates. Here, e su b = oo corresponds to a perfect 
conductor, while e su b = 7.8 and 3.12 correspond to TiC>2 and sapphire, respectively. 
The case of e su b = 1.6 is just for illustration. The interaction strength between the 
particle and the substrate is modulate by the contrast factor given by Eq. El which is 
a function of e su b- For distances, z/r< > 1, the energy goes to faster, as smaller is 
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Figure 1. Dimensionless energy as a function of z/r < for an oblate spheroid with 
r>/r< = 1.4 and different substrates. The inset in (a) is shown in (b). 



the value of e su b, as shown in Fig. ^a). At small distances, z/r < < 0.4, the interaction 
between the particle and substrate increases considerably, because the largest value of 
the multipolar interaction L becomes larger, as shown in Fig. ^b). However, the value 
of L is independent of e su b, and the observed differences at a given distance are due to 
the fact that the strength of the multipolar surface plasmon interactions are modulated 
with the factor f c , such that, for larger values of the substrate dielectric constant the 
force strength increases. In Fig. ^ we observe that the energy at a given distance is 
larger for greater values of e su b- The same general behavior of the energy as a function 
of the dielectric function of the substrate is also found for prolate spheroids. 

Now let us examine in detail the influence of the geometrical parameters on the 
energy between different spheroidal particles and a sapphire substrate, as a function 
of the dimensionless distance z/r < . In Fig. |2fa), the dimensionless energy for oblate 
spheroids with different aspect ratios between the semi- axes, r > /r < , is shown. At a fix 
distance, the energy for oblate spheroids is larger when the aspect ratio r > /r < — > 1, 
i.e., when the spheroids tend to the spherical shape. Since the substrate is always the 
same, this means that the value of L becomes larger when r > /r < — > 1 at any distance. 
In Fig. |2^b), the same is shown for prolate spheroids, however, we observe two different 
regimes in this case. When z/r < < 1.2, the value of L becomes larger when r>/r< — > 1, 
but the contrary occurs when z/r < > 1.2. This means that the value of L not only 
depends on the separation distance, but also on the axes aspect ratio r > /r < , and the 
specific symmetry of the particle. From Fig. |21 we also observe that the energy behavior 
is different for oblate and prolate spheroids, which is more evident when z/r < < 0.4. The 
energy for oblate particles becomes more negative faster than for prolate ones, because 
the value of L is larger for oblates. In conclusion, the exponent j3 of the power-law 
behavior of the energy depends on the specific geometrical parameters z, r > and r<. 

To examine in detail the energy dependence with the geometry of the system, let 




consider different spheroidal particles at a fix distance z and all with the axis normal the 
surface r±, such that, z/r± = 0.25 is a constant. In Fig. Efa) is shown the dimensionless 
energy as a function of the aspect ration r = r\\/r± between the axis parallel to the 
surface (r\\) and the one normal to it. The case when r < 1 corresponds to prolate 
particles, while r = 1 is the value for spheres, and when r > 1 we have oblate ones, 
as it is illustrated in the inset of Fig. |3Ja). Here, one can observe clearly the energy 
reliance on the geometry. For prolate spheroids the relevant multipolar interactions do 
not change dramatically with the aspect ratio of the axes. This means that the value of 
L increases smoothly as the aspect ratio goes to r = 1. On the other hand, for oblate 
particles the number of multipolar interactions is more sensitive to r, and the values of 
the largest multipolar interaction rises dramatically with small increments of r. 




Figure 3. (a) Dimensionless energy of spheroidal particles at a fix distance z, as a 
function of the aspect ratio between the axis, (b) Dimensionless energy of prolate 
spheroids at a fix z and r < constant, as a function of r> /r<. 
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In Fig. Efb), the energy of prolate particles, at a fix z and r< constant, with 
z/r > = 0.1, as a function of r > /r < is shown. In this case, the projected area of the 
particle over the substrate is always the same, since the minor axis is constant for all 
cases (see inset in Fig. Efb)). Here, we observe that when the major axis increases the 
energy decreases because the reduction of number of relevant multipolar interactions 
involved, yielding different exponent of the power law function of z for the energy. 
Therefore, we have shown that the energy is very sensitive to the geometry of the 
system. 

4. Discussion 

The above results might provide an insight into the range of validity of the Proximity 
Theorem Approximation, which was developed by Derjaguin and Abrikosova ^1] to 
estimate the force between two curved surfaces of radii R\ and i?2- The Proximity 
Theorem Approximation assumes that the force on a small area of one curved surface 
is due to locally "flat" portions on the other curved surface. Within the Proximity 
Theorem Approximation and when z <C Ri,B,2, the force per unit area between two 
curved surfaces is: 



where V(z) is the Casimir energy per unit area between parallel plates separated by a 
distance z. In the limit, when R\ — > oo and R2 = R, the problem reduces to a curved 
surface of radius R and a flat plate, yielding 



The force obtained is a power law function of z, and at "large" distances J~{z) oc z~ 3 , 
while at short distances T{z) oc z~~ 2 . The force is proportional to the radius of the 
curved surface, and to the inverse of its projected area on the flat plate. Therefore, 
particles at a given distance with the same radii of curvature would feel a force with the 
same power law function of z. Curved surfaces with the same projected area at a given 
distance would feel also a force with the same power law function of z 

In Fig.EJ we show the dimensionless energy as a function of the aspect ratio r>/r<, 
for (a) oblate and (b) prolate particles. We consider the case when all oblate (prolate) 
particles have the same curvature and are at a distance z/r > = 0.25 (0.1). The energy 
for oblate, as well as, for prolate particles is a power law function of z, which depends 
on the geometrical parameters, even when the curvature of the particle is not changed. 
In the case of prolate particles, the value L increases smoothly as the aspect ratio also 
does. On the other hand, for oblate particles the number of multipolar interactions is 
more sensitive to axes aspect ratio, and the values of the largest multipolar interaction 
increases dramatically as r>/r< also does. 
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r > /r < -V r < 

Figure 4. Dimensionless energy as a function of r>/r< for (a) oblate and (b) prolate 
spheroids with the same curvature, at a fix distance. 



5. Conclusions 

The role of geometry in dispersive forces in the non-retarded limit is studied by 
calculating the energy from the interacting surface plasmon of macroscopic bodies. We 
analyze in detail the interaction of oblate and prolate particles with a substrate. When 
the particle is close to the substrate the multipolar interactions induced by the substrate 
modify the electromagnetic response of the system. In general, we find that the energy 
is described by a power law function whose exponent depend on the minor and major 
axes of the spheroid, as well as, with the separation between bodies. 
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